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ANALYSIS OF A WATER TABLE SIMULATION OF A 
TRANS P I RAT IONALLY COOLED NOZZLE 
by Albert F. Kascak 
Lewis Research Center 

SUMMARY 

This report contains an analysis of a water table simulation of a transpirationally 
cooled nozzle. Two flow models of the water table boundary layer are formulated and 
then analyzed to predict experimental results. The approach taken is to assume that the 
pressure is constant in the transverse direction and dictated by the main flow. It is also 
assumed that the flow is inviscid and that there is no mixing between the injected flow 
and the main flow. 

The first model analyzed is one in which there is no mixing in the boundary layer, 
and the second model analyzed is one in which there is complete mixing in the boundary 
layer. The models are simple enough so that extending them to apply to a gas nozzle, 
including density differences and heat transfer, should be possible. 

The results of this analysis are presented in the form of four relations. The first 
is a relation governing the flow through any water table channel; the other three are 
relations governing the boundary layer thickness, mass, and momentum. The equations 
are coupled by the fact that the free -flow channel is equal to the physical channel minus 
the boundary layer thickness. If the flow is not choked, the height of the main flow at 
the start of injection is an independent parameter and the problem is an initial- value 
problem. If the flow is choked, the height of the main flow at the start of injection is 
dependent on the area of the throat. The area of the throat is unknown since the bound- 
ary layer is unknown. The problem then becomes a two-point boundary value problem. 

Both models predict the photographic data with similar results; but inherent two- 
dimensional and parallax errors preclude choice of the best model. The largest differ- 
ence in the two models occurs in the boundary layer thickness, with less of a difference 
in boundary layer momentum and only a slight difference in boundary layer mass. 

Transpiration flows of the order of a few percent of the main flow result in only 
slight perturbations of the main flow, while transpiration flows of the order of tens of 
percent result in major changes. In general, the boundary layer thickens rapidly near 
the start of injection because of low acceleration forces, while it thins near the throat 
because of high acceleration forces. The transpiration fluid increases the radius of 
curvature near the throat and decreases the divergence angle of the free flow downstream 
of the throat. Thus, the design of a transpirationally cooled nozzle can incorporate a 
smaller radius of curvature at the throat and a large divergence angle. 



INTRODUCTION 


Advanced propulsion devices such as the gas-core nuclear rocket have predicted 
specific -impulse capabilities between 2000 and 7000 seconds (ref. 1). Regenerative 
cooling may not be able to cool these high-specific-impulse nozzles (ref. 2). A possible 
way to cool these nozzles may be by transpiration (ref. 3). A coolant first cools the wall 
as it passes through it and then forms an insulating layer adjacent to the wall. 

The interaction of the transpiration coolant and the hot propellant is a complicated 
phenomenon. In order to get insight into this problem, a water table simulation of a 
transpirationally cooled nozzle was used for flow visualization studies by H. C. Perkins 
and R. B. Kinney of Arizona University under NASA grant NGR-03-002-213. (There is 
an analogy between compressible flow through a nozzle and flow over a water table, 
ref. 4. ) 

This study formulates two analytical flow models of the water table (with transpira- 
tional flow) and then analyzes these models to predict the experimental results. The ap- 
proach taken is to assume that the pressure is constant in the transverse direction and 
is dictated by the main flow. It is also assumed that there is no mixing between the in- 
jected flow and the main flow. The first model analyzed is one in which there is no mix- 
ing in the boundary layer, and the second model analyzed is one in which there is com- 
plete mixing in the boundary layer. 


ANALYSIS 

A water table which is used to simulate a transpirationally cooled nozzle is shown in 
figure 1. The water table is a level surface with convergent-divergent vertical side 
walls. The level surface and the side walls form a channel (open on the top) through 
which water flows. Transpirational fluid flows from a plenum outside the convergent- 
divergent channel walls, through the permeable channel walls, and into a boundary layer 
adjacent to the channel walls. All the fluid in the boundary layer is transpiration fluid 
and comes through the side channel walls. The surface separating the main flow from 
the boundary layer is assumed to be parallel to the channel wall. The height of the fluid 
in the boundary layer is equal to the height of the main stream. The flow area through 
the channel wall is limited by a liner to the height of the main stream. 

The main flow analysis is given in reference 5 and is only sketched herein to show 
the simularities between the main flow and the boundary layer. 

The flow area is 


A = bh 


(1) 


2 


(All symbols are defined in the appendix. ) For steady-state, incompressible, inviscid, 
one -dimensional flow the continuity and momentum equations are 


m = pvA = pvbh = constant 

(2) 

dM = _ £& A. (bh 2 ) + ££ h 2 ^ 

(3) 


dz 2 dz 2 dz 


Where the average pressure is 1/2 pgh and the terms on the right side of the momentum 
equation represent the forces acting. The first term represents the force on a control 
surface normal to the axis. The second term represents the axial component of the 
force on a control surface adjacent to the side wall of the channel. The control surface 
on the top of the channel is a free surface and does not contribute. The control surface 
on the bottom is level (or parallel with the z-axis) and does not contribute in the hori- 
zontal direction. Differentiating in equation (3) results in 

m — = - pghb — (4) 

dz dz 

Substituting equation (2) into equation (4) gives 

pvhb ^ = - pgbh ^ (5) 

dz dz 


Then equation (5) reduces to the standard momentum relation for a water table, 

0 = — (v 2 + 2gh) (6) 

dz 

Equation (6) can be integrated from a point where the velocity is zero (stagnation condi- 
tion) to any arbitrary point, 


v = 2g(h g - h) 


(V 


Equation (7) can be substituted into equation (2) and the results can be written as 

1 


o . (M.\* - M 


'm \ 2 


2gh 


3 \pby 


( 8 ) 
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This cubic equation has the following solutions: 


A. = i + 2 ( 9a ) 

h s 3 3 ^ 3 y 


where 


cp = cos 


-1 



27 
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(pb 


2 


(9b) 


and k is equal to 0, 1, or 2. The boundary conditions far upstream are h — h as 

o 

b — The only value of k which meets these conditions in equation (9) is k = 0. The 
boundary conditions far downstream are h — 0 + as b — 00 . The only value of k which 

meets these conditions in equation (9) is k = 2. At the throat k = 0 and k = 2 should 
give the same height. Therefore, <p must be equal to n and cos <p = - 1. Then from 
equation (9b) 


h 


s 



(10) 


Equations (7), (9a), and (10) form a complete description of a water table flow once the 
mass flow rate and channel width are specified. 

The free-flow area is shown in figure 1 by the curves. The free-flow channel width 
is approximately equal to the width of the channel walls minus the boundary layer thick- 
ness on both sides. The curved flow width is approximated by the actual channel width. 
This is a common approximation for one-dimensional nozzle flow. 


b=b w -26 (11) 

(There is no boundary layer on the bottom because there is no transpiration fluid injected 
on the bottom and the main flow is inviscid. ) It is assumed that the pressure (or height) 
is constant in the transverse direction and is dictated by the main flow (eqs. (7), (9), 
and (10)). It is also assumed that there is no mixing between the main flow and the 
boundary layer and that the boundary layer is inviscid. The continuity and momentum 
equations for the boundary layer are 


dS 


^w 

dS 


( 12 ) 
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( 13 ) 


^ = _ £g A (26h 2 ) + £& h 2 — (26) 
dS 2 dS 2 dS 


Where S is the distance measured along the channel wall and the transpiration fluid has 
no momentum in the S direction, when it initially enters the channel. The distance 
along the wall can be related to the axial distance by the geometry of the channel wall, 



The transpiration fluid can be considered a jet discharging from a higher pressure 
outside the wall to a lower pressure inside the channel wall, 



dS 


C P y2g(h p - h) h 


(15) 


where C is the discharge coefficient and hdS is the flow area. When equations (14) 
and (15) are used, the continuity and momentum equations for the boundary layer become 


■fob 

dz 


2pC V2g(h p - h) 



h 


(16) 


dMi. jl 

— - = - pgh(26) — (17) 

dz dz 


There are two models used to calculate the relation between the boundary layer 
thickness mass and momentum. The first assumes no mixing in the boundary layer; the 
second assumes complete mixing in the boundary layer. 

In the no-mixing boundary layer model; the boundary layer is composed of stream 
tubes extending from the point of injection on the channel wall and flowing adjacent to 
the wall and out the nozzle (fig. 2). Each stream tube behaves like a miniature nozzle, 
except that rather than the width being specified the pressure (or height) at each location 
is specified from the main flow calculations. Rewriting equation (8) in terms of this 
miniature nozzle and solving for its width gives 


d6 



(18) 
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The mass flow rate through this minature nozzle is equal to one -half the increment in 
the mass flow rate of the boundary layers at the point of injection, 


*,=1*6,1 = <VSi 


( 19 ) 


The subscript i represents the point of injection. The stagnation height of this minia- 
ture nozzle is just the height at the point of injection. (That is, the injected flow is nor- 
mal to the wall, and the tangential velocity is zero. The normal component of the flow 
is small since the wall area is large, and therefore any recovery of the flow in the tan- 
gential direction can be neglected. ) Symbolically, integrating equation (18) along the 
axis from the point at which injection starts Zq to the point of interest z gives 

6 = — f (^) (20 

2 J WiPh 

where the i subscripts denote the quantities changing in the integration. Likewise the 
momentum associated with this miniature nozzle is 



^b = V 2g(h s,w - h) ^w 


Applying the same reasoning and integrating give 



Equation (21) can be differentiated with respect to the axial distance and, when equa- 
tion (20) is used, results in a restatement of equation (17). Therefore, equation (21) 
is not an independent relation. A complete description of the water table with trans- 
pirational boundary layers is given by the stagnation heights upstream and outside the 
channel walls, by the geometry, and by equations (7), (9), (10), (16), (17), and (20). 

The second model analyzed is one in which there is complete mixing in the boundary 
layer. The boundary layer is characterized by velocity v^. The mass flow rate is 

= pv b (26)h (22) 


6 



and the momentum associated with the mass flow rate is 



= v b m b 


( 23 ) 


Equations (22) and (23) can be solved for the boundary layer thickness and velocity: 



2phM b 


M, 


v b = 


m 


(24) 


(25) 


Equations (7), (9), (10), (16), (17), and (24) form a complete set for this model. 

The preceding equations can be nondimens ionalized if a* is defined to be a refer- 
ence length and /3 is a constant equal to 

/3 ^-i= ( 26a ) 

2pa*h s y2gh s 


The parameters are then defined by 


e 


y = 


h6 


/3h g a* 


4 = 


V 


w 


/3a* 

b w 
2 /3a* 

m. 


m 


(26b) 

(26c) 

(26d) 

(26e) 

(26f ) 
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u = 


(26g) 


M b 

y2gh s m 

1 = — 

/ 3a * 

Substituting equations (11) and (26) into (8) yields 

e =y v - — L= 

V 1 - y 

Substituting equations (26) into (16) and (17) yields 


(26h) 


(27) 


~ = Cy y Yp - y ^1 + ( r /) 2 (28) 

^ = - £ & (29) 

d£ 2 d£ 

where tj the channel geometry, C the discharge coefficient, y p the stagnation height 
outside the channel wall are given functions of £. Substituting equations (26) into equa- 
tion (20) yields the nondimensional relation for the no-mixing boundary layer model, 



Substituting equations (26) into equation (24) yields the corresponding relation for the 
completely mixed boundary layer model. 


e 



u 


(30b) 


Equations (27) to (30) form a complete set of four equations in four unknowns. The 
initial conditions for this set of equations are that, at the start of injection of the trans- 
pirational fluid into the boundary layer, the boundary layer thickness, mass, and mo- 
mentum are zero: 


8 


<(« 0 > = ° 


(31a) 


wu 0 ) = o 

(31b) 

uU 0 ) = 0 

(31c) 


If the flow is not choked flow, the height of the main flow at the start of injection is an 
independent parameter 


yU 0 ) = y 0 (3 id) 

and the problem is an initial-value problem. If the flow is choked, the height of the main 
flow at the start of injection is dependent on the area of the throat. The area of the 
throat is unknown since the boundary layer is unknown. The problem is a two-point 
boundary value problem. 

In order to avoid the singularity in equation (30a) and to have finite derivatives at 
the start of injection, the independent variable is changed from £ to y. (This assumes 
that y is a monotonic function of |.) Therefore, £, e, w, and u are functions of y. 
Equations (28), (29), and (30a) become 


W* = Cy yy p - y V 



The initial condition (3 Id) comes from equation (27) when e equals zero 


vU 0 ) = 


1 

yo V 1 - *o 


(32) 

(33) 


(34) 


(35) 


Assume that near the start of injection 

e ~ (y 0 - y) a 
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where a is greater than zero since eg is equal to zero. From integrating equa- 
tion (33) 


u ~ (y 0 - y ) 0+1 

When equation (30b) or Abel's solution of equation (37) is used, 

W ~ (y Q - y) Q ' +(1/2) 

From integrating equation (32) 


4 - 4 0 ~ (y 0 - y) a ' +(1/2 ' 

Differentiating equation (27) and substituting yields a equal to 1 as the only finite solu- 
tion. Therefore, £q, Wq, and Ug are zero and 


£ , 2 -% 
2y 0 ^-y 0 ) 3/2 

Equation (36) was derived by differentiating equation (27) 


(36) 


e' = 7j + y n't' -- (i - y)" 3/2 
2 


(37) 


and using the initial values. 

Equation (34) is a form of Abel's Integral Equation; and as such, it has an inverse 



Numerically, it is convenient if an additional independent variable transformation is 
made, 




( 39 ) 


Then equation (38) becomes 
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w* = + - 


The integral in equation (40a) exists if e' is finite. Equation (30b) can be differentiated 
to give the completely mixed model equivalent of equation (40a), 


W‘ = 


Equations (28), (29), (37), and (40) form the equation set which is solved numerically. 
The analytical solution near the start of injection is 


e = eo ( y ' y 0> 


u = - — (y - y 0 ) 2 

4 u 


t - r-JH r -^ 2 % - y > 3/2 

c o y oV y p,o- y o V lt( w 

w = - B «o(y 0 ‘ y)3/2 


where B for the no-mixing case is 


B = — 


and B for the complete -mixing case is 


B = - 


Equation (40) can be written as 


W' = Fe 1 + G 
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where F and G are auxiliary functions. For the completely mixed boundary layer 
model the auxiliary functions are 


F 


u 

2W 


and 


G = - 



(43a) 


(43b) 


For the no-mixing boundary layer model and when the trapezoidal method is used to ap- 
proximate the integral, the auxiliary functions become 



(44a) 


and 


i / 2 i 


2 vy 0 -y 

£ » dXi+£ ’ 


y-Ay 


(44b) 


Equations (32), (37), and (40) can be solved simultaneously for and the result is 


r = 


2Ft? - F(1 - y)~ 3//2 + 


2G 


2C y yy - y ^1 + ( V ') 2 - 2Fy77* 


(45) 


By using the initial conditions (eqs. (36) and (41)) the differential equations (32), (33), 
(37), and (45) can be solved numerically by the "Predictor -Corrector Method. " 


DISCUSSION 

The analysis of a water table simulation of a transpirationally cooled nozzle is given 
in the previous section of this report. The results of this analysis are in the form of 
four relations. The first is a relation governing the flow through any water table chan- 
nel, and the next three are relations governing the boundary layer thickness, mass, and 
momentum. The equations are coupled by the fact that the free-flow channel is equal to 
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the physical channel minus the boundary layer thickness and also by the fact that the 
pressure (or height) in the boundary layer is the same as in the main flow. Two bound- 
ary layer models were assumed; first, no mixing and, second, complete mixing in the 
boundary layer. 

In order to check the validity of this analysis, a calculation was made for the water 
table of Perkins and Kinney. A "plan view" of the nozzle walls for this water table is 
shown in figure 3. The main flow and the flow through sections 1, 2, and 3 were me- 
tered. The flow through sections 1 to 3 was dyed and, when photographed, appeared 
darker than the undyed main flow. Experimental results were recorded photographically 
for the case of the main flow being equal to 2. 27 kilograms per second (5. 0 lbm/sec) and 
the flow through sections 1, 2, and 3 being equal to 0. 227, 0. 227, and 0 kilograms per 
second (0. 5, 0.5, 0 lbm/sec), respectively. 

For this case the calculation was done in the following fashion: First, it was as- 
sumed that there was no boundary layer adjacent to the channel wall, (e =0). The height 
of the flow y was calculated by using equation (27). With this height, the mass flow 
rate in the boundary layer W was calculated by using equation (32). These two results 
were then substituted into equation (34) and a new boundary layer thickness e was cal- 
culated. The process was repeated, with the new value of the boundary layer thickness, 
until there was no change from one iteration to the next. An instability arose in the vi- 
cinity of the throat and a smoothing technique was applied so that the position of the 
throat (critical flow) was stable. 

Figure 4 shows the plan view of the water table nozzle. It shows the position of the 
channel wall and the adjacent boundary layer as a function of axial position, from the 
start of injection down past the physical throat. The flow is choked, and as such the 
physical throat and critical flow do not necessarily occur at the same axial position. In 
general, the boundary layer thickens rapidly near the start of injection and thins some- 
what in the high-acceleration region near the throat. The transpiration fluid in general 
increases the radius of curvature near the throat and decreases the divergence angle 
downstream of the throat. This suggests that the physical radius of curvature at the 
throat can be made much smaller and the physical divergence angle can be made larger 
for a transpirationally cooled nozzle than for a conventional nozzle. 

Figures 4(a) to (f) show the results for various iterations. The change is small for 
the first two iterations and is negligible for the next four iterations. On the sixth itera- 
tion, instabilities build up (even though a smoothing routine is used) and cause numer- 
ical difficulties in the vicinity of the throat. 

These calculations were compared to the experimental results shown in figure 5. In 
the photograph there is some problem with parallax, and therefore only general trends 
can be considered. In general, the calculated boundary layer is in close agreement with 
the photographed boundary layer. The region of largest disagreement is just upstream 
of the throat and is probably caused by two-dimensional effects. 
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The preceding calculations were for the no-mixing boundary layer model. Similar 
calculations were made for the completely mixed boundary layer model, and these re- 
sults are very close to those shown in figure 5. Figure 6 shows the difference between 
the two models near the start of injection. The boundary layer thickness, mass, and 
momentum are plotted as a function of axial distance from the start of injection. In gen- 
eral, both models give similar results. The largest difference occurs in the boundary 
layer thickness. A smaller difference occurs in the boundary layer momentum; and no 
difference occurs in the boundary layer mass. These results are, in general, also true 
farther downstream, as shown in figure 7. 

Figure 7 shows the results of the analysis for a simplified geometry. The geometry 
considered was straight convergent channel walls at a angle of 45° to the centerline. 

The channel walls, if extended, would have met at the origin of the axial coordinate. 

The channel wall discharge coefficient was assumed to be 0. 01. The stagnation height 
outside the channel wall was assumed to be equal to the main-stream stagnation height. 
The problem was an initial-value problem in which the initial nondimensional height was 
assumed to be 0. 95, 0.90, 0.80, and 0.70. The problem was solved from the initial 
point to the point of choked flow for both models considered in the analysis. 

Figure 7(a) shows the nondimensional mainstream height as a function of nondimen- 
sional axial distance. The main-stream height is essentially the same for the four dif- 
ferent starting conditions. In this case the boundary layer does not influence the main 
stream very much, or the flow is weakly coupled. The reason for this is that the bound- 
ary layer thickness is small compared to the channel width. This is shown in figure 7(b). 

Figure 7(b) shows the nondimensional boundary layer thickness as a function of non- 
dimensional axial distance. From the assumed geometry the channel width is equal to 
twice the absolute value of the axial distance. The largest boundary layer is about 
1 percent of the channel width, and thus has a small effect on the free-flow area. The 
boundary layer thickens rapidly near the start of injection and thins rapidly near the 
throat. The reason for this is shown in figure 7(c). 

Figure 7(c) shows the nondimensional boundary layer mass as a function of nondi- 
mensional axial distance. The curves have almost a constant slope and therefore almost 
constant mass added to the boundary layer per unit axial length. Near the start of injec- 
tion, the acceleration (the slope of the height-against-distance curve) is low, and there- 
fore the boundary layer must thicken to accommodate the additional mass. Near the 
throat the acceleration is large, and therefore the boundary layer thins. This accelera- 
tion near the throat is shown in figure 7(d). 

Figure 7(d) shows the nondimensional boundary layer momentum as a function of 
nondimensional axial distance. Near the start of injection the change in momentum is 
small. Near the throat the change in momentum is large. 

Figure 8 shows a strongly coupled case. That is, the boundary layer greatly affects 
the main flow. The case considered was the same as in figure 7, except that the dis- 
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charge coefficient was increased from 0. 01 to 0. 1. This resulted in more mass in the 
boundary layer and therefore a thicker boundary layer (which in turn decreased the free- 
flow area). 

In summary, this report presents two simple models for a water table nozzle, in- 
cluding transpirational flow. The models should be applicable in a region where pres- 
sure forces dominate over viscous forces (near the throat). The models are simple 
enough so that extending them to apply to a gas nozzle, including heat transfer and den- 
sity differences, should be possible. 


CONCLUSIONS 

Analysis of a water table simulation of a transpiration-cooled nozzle is given in this 
report. Two boundary layer models were assumed; first, no mixing and, second, com- 
plete mixing in the boundary layer. The analysis was compared to photographic data and 
resulted in the following conclusions; 

1. The two boundary layer models predicted the data within 20 percent over most of 
the range. 

2. Both models predicted similar results; but inherent errors precluded a choice as 
to which model fit the data best. 

3. The largest difference between the data and the analysis were caused by two- 
dimensional effects being omitted from the analysis and by parallax in the photographic 
data. 

4. The largest differences in the two models occurred in the boundary layer thick- 
ness. Less differences occurred in the boundary layer momentum; and only a slight 
difference occurred in boundary layer mass. 

5. Transpiration flows of the order of a few percent of the main flow resulted in 
only slight perturbation of the main flow. Transpiration flows of the order of tens of 
percent of the main flow resulted in major changes of the main flow. 

In general, the transpiration fluid increased the radius of curvature near the throat 
and decreased the divergence angle of the free flow downstream of the throat. Thus, the 
design of a transpirationally cooled nozzle can incorporate a smaller radius of curvature 
at the throat and a larger divergence angle. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 26, 1973, 

503-04. 
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APPENDIX - SYMBOLS 


o 

A flow cross-sectional area, m 

a* reference length, m 

B constant defined in eq. (41e) or (41f) 

2 

b width of main flow, m 

C channel wall discharge coefficient a function of axial position 
d nondimens ional boundary layer thickness 
F auxiliary function defined in eq. (43) or (44) 

G auxiliary function defined in eq. (43) or (44) 

O 

g acceleration of gravity, 9.8 m/sec 
h height of flow, m 
k constant equal to 0, 1, or 2 

2 

1)1 momentum flow rate, kg-m/sec 
m mass flow rate, kg/sec 

5 curved distance along channel wall, m 
u nondimens ional momentum flow rate 

v velocity, m/sec 
W nondimens ional mass flow rate 
x auxiliary independent variable 
y nondimens ional flow height 
z axial distance, m 
ot exponent used in initial conditions 
nondimensional constant 

6 boundary layer thickness, m 

e nondimensional boundary layer flow area 
7 ] nondimensional channel wall width 
£ nondimensional axial distance 

Q 

p density, kg/m 

cp angle used in solution of cubic equation 
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Subscripts: 

b boundary layer 

i injection point 

p plenum 

s stagnation condition 

t throat condition, choked flow 

w wall 

1, 2, 3 sections 1, 2, 3 
0 start of injection 

Superscripts: 

(') differentiation with respect to given independent variable 
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(e) Iteration 5. (f) Iteration 6. 

Figure 4. - Plan view of water table nozzle, showing position of channel walls and adjacent boundary layer. Flow rate, kg/sec (ibm/sec): 
main, 2. 27 (5); section 1, 0. 227 (0. 5); section 2, 0. 227 (0. 5); section 3, 0 (0). 










Figure 6. - Boundary layer thickness, mass, and momentum 
as function of axial distance from start of injection of 
transpiration flow. 
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